Control of laser wake field acceleration by plasma density profile 
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We show that both the maximum energy gain and the accelerated beam quality can be efficiently 
controlled by the plasma density profile. Choosing a proper density gradient one can uplift the 
dephasing limitation. When a periodic wake field is exploited, the phase synchronism between the 
bunch of relativistic particles and the plasma wave can be maintained over extended distances due 
to the plasma density gradient. Putting electrons into the n— th wake period behind the driving 
laser pulse, the maximum energy gain is increased by the factor 2im over that in the case of uniform 
plasma. The acceleration is limited then by laser depletion rather than by dephasing. Further, we 
show that the natural energy spread of the particle bunch acquired at the acceleration stage can be 
effectively removed by a matched deceleration stage, where a larger plasma density is used. 

PACS numbers: 52.38.Kd,52.65.Rr,52.27.Ny 



Plasma-based schemes of electron acceleration have re- 
cently demonstrated impressive progress. Quasimonoen- 
ergetic electron bunches with the energy up to 1 GeV and 
with the charge of 50 pC have been generated in exper- 
iments All-optical methods for control of the bunch 
parameters have been developed also [2(. It is generally 
believed that electrons in these experiments have been ac- 
celerated in the "Bubble regime" [3(. At the same time, a 
significant advance in plasma profile engineering has been 
achieved to make laser plasma interaction more efficient: 
plasma capillaries for laser guiding 0] , fabrication of cor- 
rugated plasma structures |f| , plasma machining Q . 

One of the main limitations on energy gain in laser- 
plasma accelerators comes from the dephasing. The ve- 
locity of relativistic electrons is slightly higher than the 
phase velocity of the wake, which is determined by the 
group velocity of the driving laser pulse. The accelerated 
electrons slowly outrun the plasma wave and leave the 
accelerating phase. 

The limitation caused by the dephasing can be over- 
come by employing a proper plasma gradient The 
profile of the plasma density should be such that the ad- 
vance of the accelerated electrons matches the advance of 
the plasma wave. The equation for plasma density profile 
in the ID configuration is 
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where oj^ix) = Aite 2 n{x) /m is the squared plasma fre- 
quency, n(x) is the plasma density, $„ = const is 
the phase of ultrarelativistic electrons trapped n plasma 
wavelengths (A p = 2irc/uj p ) behind the laser pulse, v gr 
is the group velocity of the laser pulse, c is the speed 
of light, e and m are the electron charge and mass 
respectively. For weakly relativistic laser pulses with 



eA/mc 2 -C 1 and for rarefied plasmas n/n c -C 1 
•g r /c ~ 1 — n(x)/2n c and 7 7 gr , where 



we can assume v 



7 is the relativistic gamma-factor of the accelerated elec- 
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mw 2 /47re 2 is the critical 
plasma density and u> is the laser frequency. The solution 
of Eq. ([T]) for the phase synchronism in laser wake field 
acceleration is 
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where tiq = n(x = 0). It follows from Eq. Q that 
the plasma density increases along the pulse propagation 
and the acceleration distance is limited by L m d since the 
plasma density goes to infinity at x — Li n h- 

Consider a short circularly polarized laser pulse with 
the Gaussian envelope a 2 (£) = a 2 , exp (— r 2 /T 2 ), where 
t = t — f dx' /vg r (x'). For simplicity, the ion dynam- 
ics, the thermal motion and the transverse dynamics of 
plasma electrons are neglected. It is also assumed that 
the accelerated electrons do not affect plasma wake struc- 
ture. The accelerating force on the relativistic electrons 
in the plasma wake can be presented in the form Q 
F x = - {^/2)a 2 ,mcLu 2 (x)Texp [-w 2 (x)T 2 /4] cos$(ic) 

, where = w p (x)r(x) ~ \Jn{£)/n c J n(£)/ (2n c ) d£ 
and £ = lox/c. For cos$ = —1 the accel- 
erating force achieves a maximum value of F x ~ 
mcujp^it /2a 2 , exp(— 1/2) when T = \f2/w p . Therefore 
for the given duration of the laser pulse there is the opti- 
mal density (uj 2 T 2 = 2), at which the accelerating force 
peaks. When the laser pulse propagates in inhomoge- 
neous plasma the pulse duration will be soon out of op- 
timal value. In addition, the nonlinear dynamics of the 
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FIG. 1: Electron acceleration in plasma layers in the frame 
of light speed (schematically). The ultra-relativistic electron 
bunch (black rectangle) is always located at the peak value of 
the accelerating field. 



laser pulse during propagation can significantly modify 
the pulse envelope and carrier frequency. The both ef- 
fects influence the energy gain and should be taken into 
account to find the plasma density profile for optimal 
acceleration. 

The electron energy gain is AS — (J F x dx). Intro- 
ducing the parameter /i = J (n/n c ) dx the phase can be 
rewritten in the form $ = and the energy gain 

can be presented as follows 



F(£,fJ,,fi) = ^-ag/t/3(C) exp 
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where the evolution of the pulse duration is assumed to 
be known /?(£) = u)T(£). Considering AS \jj>(£)] as a 
functional the Euler-Lagrange equation can be derived 
for /i (£) providing the peak gain 
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Assuming that the laser pulse is short, uj p T <C 1, so that 
the pulse duration effects can be neglected /?(£) = 0, 
Eq. (J6J) reduces to Eq. (P) with <!>„ = — tt(1 + 2n), where 
n = 1,2,... 

Integrating Eq. ([1]) for homogeneous plasma (n = tiq — 
const) the wake phase at the bunch position can be cal- 
culated <I> = -37r/2 + (n /n c ) 3 / 2 £/2, where it is as- 
sumed that = 0) = — 37r/2. The electron accel- 
eration is possible in the range < x < £hom, wnen 
-7r(3/2 + 2n) < $ < -7r(l/2+2n), where n = 1,2,... and 
Lhom = 2TTc(no/n c )~ 3 / 2 /uj is the well known expression 
for detuning length in homogeneous plasma Q . Integrat- 
ing Eq. Q we obtain the known expression for the peak 
gain in the electron energy for homogeneous plasma Q 
Afhom = 2\/27rexp(— l/2)a§(n c /no)mc 2 , where optimal 
duration of the laser pulse T = ^/2/ui p was assumed. 

For a short laser pulse u p T -C 1 the peak gain can be 
achieved when the wake phase at the electron position is 
<!>„ = — 7r(l + 2n) and the plasma density profile obeys 
Eq. |2]). Integrating Eq. ([4]) the energy gain takes the 
form 



AS m h(x) ~ So 
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where So = (i/2)^a^mc 2 (no/n c )(ujT)(uiLi n h/c). The 
energy gain over the distance < x < Li n h is 
ASinh ^ S o ~ A£hom2~ 3/2 cxp(l/2)w p oT<I>„, where 
uip = 4ire 2 no/m. It follows from the obtained expression 
that the gain increases as $„. Therefore, the electron ac- 
celeration is most efficient when the electrons are loaded 
at the peak accelerating field as far behind the laser pulse 
as possible. For arbitrary values of to p T, and the phase 
synchronism $„ = const ensured by the plasma profile 
dU), the energy gain is 



AS inh (x) = So [*(x) - *(0)] , 
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where S(x) = uo 2 p (x)T 2 /A and 6(0) = u 2 pQ T 2 /<i = 
n uJ 2 T 2 /4n c . In the limit lo p T « 1 (5 « 1) Eq. © 
is reduced to Eq. ||7J). 

Keeping the bunch always in the same wave period 
behind the laser pulse may become unpractical, because 
the plasma density would vary too strongly over the full 
acceleration stage. To optimize the process we propose a 
layered plasma density profile. The acceleration scheme 
is illustrated in Fig. [T] Let the electrons be trapped in 
the n-th plasma period during acceleration in the first 
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FIG. 2: Electron bunch acceleration and energy spread re- 
duction in plasma layers: (a) the plasma density profile, (b) 
the mean energy of the electron bunch (c) electron energy vs 
electron position in the bunch. Electrons are accelerated by 
the first laser pulse in the first two layers (0 < x < 43370 Al) 
whereas they are decelerated in the third and fourth layers 
(43370A L < x < 49820Al). The solid line and the dotted line 
correspond to the PIC simulation results and the theoretical 
estimates, respectively. To compare, the dashed line shows 
acceleration in homogeneous plasma with the constant den- 
sity no = 0.001n c . The energy distributions in frame (c) are 
shown at the beginning of the deceleration at x = 43370 Xl 
(line 1), at x = 45370 X L (line 2), at x = 47370 X L (line 3) 
and the end of deceleration at x — 49820 Xl (line 4), respec- 
tively. The laser pulses are circularly polarized with Gaussian 
profile and wavelength 1 /im. The laser pulse parameters are 
T — 10.6A_l/c, ao = 0.6 for the pulse on the acceleration stage 
and T = 5.3Xl/c, ao = 0.5 for the pulse on the deceleration 
stage. 



The density value of n 2 = n/$^_ x /$^ is chosen such that 
the electrons in the second layer become located in the 
(n— l)-th plasma period behind the laser pulse. The wake 
phase at the electron bunch center is = w(2n — 3) 

that again corresponds to the maximum of the accelerat- 
ing force. Inductively, the plasma density increases from 
n m = n f&n-m/ ^n-m-i t° n f m the m-th plasma layer 
in accordance with Eq. ([2]) and the electrons are trapped 
in the (n — m)-th plasma period behind the laser pulse. 
The successive acceleration in the several plasma layers 
leads to a significant increase in the electron energy gain 
and can be accomplished by the same driving laser pulse 
until it is depleted. The transition between the plasma 
density layers can be smooth, but much shorter than the 
length of the layers themselves. 

Further, the proposed scheme allows us to control the 
energy spectrum of the bunch electrons. Naturally, elec- 
trons located at the bunch center get the greatest energy 
gain while electrons located at the bunch edges get the 
smallest energy gain. This leads to an energy spread of 
the bunch after the acceleration stage. To reduce this 
spread we suggest to use a matched decelerating stage 
with a higher plasma density. Although the decelerat- 
ing stage will reduce the net energy gain, its effect on 
the energy spread reduction is much stronger, because 
the plasma wake period decreases as the plasma density 
increases. 

To get an estimate of the energy spread reduction we 
assume that the electron bunch duration is much shorter 
than the plasma wave period while the plasma density 
profile and energy gain are described by Eqs. and (O, 
respectively. It is also assumed that the bunch center is 
always synchronized to the maximum of accelerating or 
decelerating forces $„ = — 7m, n = 1,2,... Expressing 
the phase $„ in uj p rj/c <C 1, where rj is the distance 
from the electron to the bunch center, the force acting 
at the electron located at the position r\ can be Taylor 
expanded as F x = F x q(1 — (w p 7//c) 2 /2), where F x q is the 
force on the electron in the bunch center. The difference 
in the accelerating force acting on electrons located at 
different positions in the bunch leads to the energy spread 
of the bunch electrons. Using Eq. ([7|) the energy gain as 
a function of the electron position in the bunch is 



plasma layer where the plasma density increases from 
n\ to some final value nt in accordance with Eq. f2j). 
The wake phase is synchronized to the electron bunch 
in this plasma layer due to the plasma density gradi- 
ent. The wake phase at the electron bunch center is 
= 7r(2rt — 1) that corresponds to the maximum of the 
accelerating force. When the plasma density reaches the 
value of n/ , the electron bunch enters the second plasma 
layer where the plasma density increases from ni to n/. 



A£(x,rj) = A£ mh (x) - a (x) r\ 
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where A£i n h is defined by Eq. (|7|). Introducing the 

energy spread as a — [(£ 2 ) — (£) ) , where (£) = 

Jj/a (£ 2 ) = Jl b b ^{£ 2 /h)d v and l b /c is the 

bunch duration, we find the evolution of the bunch en- 
ergy spread 
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Let the bunch be accelerated in the first plasma layer 
where the plasma density increases from n-y to Any. Then, 
it is decelerated in the second layer where the density in- 
creases from rii to 4n2- The parameters of the laser pulse 
are assumed to be the same in the both layers. The to- 
tal energy gain after passing the two plasma layers is 
S = £i[l - ($2/ < i ) i)(^i/^2) 1/2 ] while the total spread 
is a = <7i|l — ( < i ) 2/ < i > i)(n2/ni) 1 / 2 |, where £\ and o\ are 
the energy gain and spread after passing the first layer, 
respectively, $1 and $2 are the phase of the bunch cen- 
ter in the first and second layers, respectively. Choosing 
nil rix = ®i/&2 the final spread is removed completely in 
this approximation. At the same time, the energy gain 
is S = £i(l — ni/ni). For n^/rii ~ 4 the bunch loses 
only one quarter of its energy after passing the deceler- 
ating layer while the energy spread will be completely 
removed. 

In order to check the validity of our simplified model, 
we have performed ID PIC simulations, using the code 
Virtual Laser Plasma Laboratory [lfj . The code has been 
supplemented with adaptive scheme: the plasma density 
is varied in each time step so that the bunch center is 
always located in the maximal value of the accelerating 
force. The plasma wake excited by the first laser pulse 
in the first two layers accelerates the electron bunches 
while the bunch is decelerated in the second two layers 
where the plasma wake excited by the second pulse (see 
Fig. [2k)- The laser pulses are circularly polarized with 
Gaussian profile. The laser wavelength Al — 1 /xm. 

The duration of the first pulse is T = 10.6Al/c and 
ao = 0.6 while the duration of the second pulse is 
T = 5.3Al/c and ao = 0.5. The plasma density at the 
beginning of acceleration is no = 0.001 n c . The bunch 
center in the first layer is located at 3.5 plasma peri- 
ods behind the pulse center (phase of bunch center is 
<I>i = — 77r). The second layer is chosen such that the 
bunch is located at 2.5 plasma periods behind the pulse 
center ($2 = — 57r) there. For the decelerating layers 
$3 = — 8tt and $4 = —6w in the third and fourth layers. 
The electron bunch was initially monoenergetic with en- 
ergy 50 MeV and duration T), = 2Al/c. 

The bunch energy achieves about 1.4 GeV with the en- 
ergy spread of about 2% after acceleration in the first two 
plasma layers (see Fig. [2j3,c). On the deceleration stage, 
the bunch energy reduces to 1.2 GeV whereas the energy 
spread reduces to less than 0.5% (see Fig. [2b, c). The 
total distance of the bunch acceleration and conditioning 
is about 5 cm. It is seen from Fig. [2^,,b that the plasma 
density profile and the bunch energy evolution obtained 
in PIC simulation agree fairly good with the theoretical 



predictions ([2]) and (JSJ. To achieve a better agreement 
we take into account the pulse compression in the second 
half of the first layer (we suppose that a = 0.7 and T — 
7.8Al/c), in the second layer (ao = 0.8 and T — 5.9Ai/c) 
and in the fourth layer (ao = 0.6 and T = 3.7Al/c) in ac- 
cordance with PIC simulation results. The deviation of 
the numerical results from theoretical estimates is caused 
by the complex nonlinear dynamics of the laser pulse dur- 
ing propagation. According to Fig. [2Jd, the energy gain 
in the inhomogeneous plasma is about 5 times more than 
that in the homogeneous plasma with n = no = const. 

The approach reported above is one-dimensional and 
thus does not take into account transverse dynamics of 
the laser pulse and electrons. The laser pulse can be ef- 
ficiently guided by plasma channels over many Rayleigh 
lengths 0, [ll[ . More complicated can be the transverse 
dynamics of accelerated electrons. It is well known Q 
that the simultaneous accelerating and focusing of elec- 
trons occurs over a quarter of the plasma wavelength. 
Another accelerating quarter of the plasma wave defo- 
cuses electrons. At the maximum of the accelerating 
fields, however, the transverse fields vanish. The use of 
a preformed plasma channel can significantly extend the 
region where the accelerating and focusing phases overlap 

In conclusion, we have proposed to control both the 
maximum energy gain and the accelerated beam quality 
by the plasma density profile. Choosing a proper density 
gradient one can uplift the dephasing limitation. The 
natural energy spread of the particle bunch acquired at 
the acceleration stage can be effectively removed by a 
matched deceleration stage, where a larger plasma den- 
sity is used. 
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